“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1966 


Verification of Mazeika's method of 
thermocline depth prediction for the 
Northeast Pacific Ocean 


Ciboci, John William 


http://ndl.handle.net/10945/9670 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 


1966 
_ CIBOCI, J. 





VERIFICATION OF MAZEIKA'S METHOD OF 
THERMOCLINE DEPTH PREDICTION FOR THE 
NORTHEAST PACIFIC OCEAN 


JOHN WILLIAM CIBOCI 


» L-BRARY 
JAVAL FOSTGRADUATE SCHOGL 


JONTEREY, CALIF. 93940 


a, KNOX LIBRARY 
WAVAL POSTGRADUATE SCHOO! 
MONTEREY CA 93¢53-510} 








This document has been approved for public 
release and sale; its distribution is unlimited. 














VERIFICATION OF MAZEIKA'S METHOD 
OF THERMOCLINE DEPTH PREDICTION 
FOR THE NORTHEAST 
PACIFIC OCEAN 
by 
John = espe eae 


Lieutenant, United/States Navy 
B. S., University of Idaho, 1959 


Submitted in partial fulfillment 
for the degree of 


MASTER OF SCIENCE IN OCEANOGRAPHY 
from the 


UNITED STATES NAVAL POSTGRADUATE SCHOOL 
October 1966 


: 
ABSTRACT 

Mazeika's method for forecasting mixed-layer (thermocline) depth 
of the upper ocean layers is discussed along with a newer version of 
this method developed by James. Using Mazeika's method primarily, a 
verification for the Northeast Pacific Ocean was completed with data 
from Ocean Weather Stations PAPA (50N, 145W) and NOVEMBER (30N, 140W) 
and a point named MIDPOINT (40N, 140W). The results indicate Mazeika's 
method is successful at Station PAPA more than seventy-five percent of 
the time during the heating season followed by a rapid decline as the 


cooling season begins. The method should be useful in the entire 


Central Subarctic Domain as described by John P. Tully. The method fails 


at NOVEMBER and MIDPOINT producing less than thirty percent success in 
prediction. James’ version did not improve the results obtained at 
Station NOVEMBER. This failure appears to be due to the controlling 
parameters for processes in the Subtropic or Transitional oceanographic 
regions (which include NOVEMBER and MIDPOINT); these differ from para- 
meters controlling oceanic processes in the Pacific Subarctic region 
(Station PAPA), which resemble those involved in the Atlantic region 
for which Mazeika's method was developed. 

Climatology data which can be used to obtain surface and 400-foot 
level temperature are also tested. The results indicate these data are 
very useful and accurate in determining the stability index required of 


Mazeika's method. 
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Le Background. 

The United States Navy is expending much effort in the direction 
of anti-submarine warfare. The proper use of anti-submarine detection 
devices hinges on the capability of correctly predicting the thermo- 
cline depth in the various oceans. Any reliable method which can be 
easily and properly used by the men and equipment of the anti-submarine 
forces would add to the overall defense capability of our country. To 
this end, the testing of one particular method of thermocline depth 


prediction was undertaken. 


a Introduction. 

Mazeika's method utilizes a subjective but realistic approach to 
the thermocline depth prediction problem. The thermocline depth pre- 
dicted by this procedure is the mean depth of the interface between 
the mixed layer and the thermocline. The method requires the. use of 
easily obtainable inputs such as forecast wind speed, fetch, wind 
duration and the synoptic weather picture. Subjective aspects of this 
method come to light in certain choices which are made by the user of 
the prediction technique, as will be explained later. 

Various types of vertical mixing are included in Mazeika's method. 
They are mechanical mixing, which is due to the combined effect of wind 
waves and associated wind currents; instability mixing due to advection 
and density changes; and mixing due to vertical convergence-divergence 
caused by changes in the overlying atmospheric pressure field, or by 
convergence-divergence in the surface current field. Mechanical mix- 
ing is introduced into the method by use of a sea state parameter which 
is based on the wind wave characteristics of the area under investiga- 
tion. The use of a stability index and corrections for the salinity 
gradient were chosen to satisfy the requirements of instability mix- 
ing. Finally, convergence-divergence is accounted for in the construc- 
tion of graphs used to determine the mixed~layer depth. 

Mazeika combined the use of these mixing parameters with the study 
of BT data, weather charts and simultaneous weather and wave observa- 
tions from the Atlantic Ocean Weather Station (OWS) CHARLIE (52N, 35W) 
to develop his prediction procedure empirically. After development was 
completed, data from Atlantic OWS BRAVO (56.5N, 51W), ECHO (35N, 48W), 
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and DELTA (44N, 41W) were used to test the method. No data from the 
Pacific ocean weather stations were used for the empirical develop- 
ment of the method or for its testing. 

Mazeika has stated that his prediction method can generally be 
used anywhere in the ocean without previous information on the mixed 
layer and thermocline; the only requirement is knowledge of the surface 
wind waves and the stability index. The testing of Mazeika's method in 
the Pacific Ocean would serve to substantiate this claim by Mazeika or 


possibly bring to light weaknesses in the method which may be present. 
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3. Mazeika's Method. 

Mazeika (1960) considered the mixed-layer thickness, h, to be a 
function of wave amplitude A, wavelength A, and a mixing -parameter k, 
which is valid at the interface between the mixed layer and the thermo- 
cline. The mixing parameter cannot be expected to remain constant. 
Since it applies to the bottom of the mixed layer, it must be dependent 
on mixed-layer thickness and on the wave parameters. [In Mie ok, its 
value depends on the stability in the thermocline. 


The relation describing orbital motion in trochoid waves 
21r 
k = Ae Fe (1) 


_includes most of the parameters, if h is considered to be a fixed value 
under given surface conditions and if stability in the thermocline is 
constant. The above function, k(A,h, A), was adopted by Mazeika in 
view of the generally accepted theory that orbital motion of particles 
due to waves decreases exponentially with depth and is proportional to 
the wave parameters. 


Solving (1) for h, 


h= Gln (R) “ 


Equation (2) can be used to compute h, if the wave parameters and 
the corresponding value of k for existing stability in the thermocline 
are known. 

In the above equations the amplitude (A) is one-half of the sig- 
nificant wave height (Hy 73)3 wave length (/\) was computed by the re- 


lation A= 3.4177 


_— where T ax is the period of maximum energy of the 


spectrum for a fully-developed sea. 
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On the assumption that vertical temperature difference is propor- 
tional to the mean temperature gradient in the upper part of the thermo- 
cline, Mazeika adopted the difference in the entire thermocline as the 
index of stability (a salinity correction could be added, see Table I). 
Since most BT data usually reach a depth of 400 feet, Mazeika selected 
the difference between the surface and the 400 foot level as the stabil- 


ity index 


OT = T, - Ty00° 


A deeper layer (e.g., 600 feet) could be taken if data were available. 

A certain initial stability exists with a given AMT in the absence 
of a mixed layer. If the mean stability index remains constant, stabil- 
ity increases from the initial value with increasing mixed-layer thick- 
ness. Increase of stability is slow when the mixed layer begins to 
form and gradually becomes greater as the mixed-layer thickness and AT 
increase. 

To connect mixing parameter values with stability and wave para- 


meters, another parameter (sea state parameter) was introduced by Mazeika 


We? n /geos 


T\_is used as a measure of sea state for wind waves. k values deter- 
mined from equation(1)with known wave parameters and mean mixed-layer 
thickness (h) values were plotted against TY values computed with 
corresponding wave parameters for eight groups of AT. Curves which fit 
the distribution of points for each AT group were computed and resulted 


in the solid-line curves of figures 1-8. h values which were related 


2 


to the k(T\) curves were then computed and plotted as the dashed-line 
curves on figures 1-8. The central h curve applies to normal mixing, 
the upper h curve applies to convergent situations with horizontal 
convergence occurring in the upper ocean layers adjacent to the point 
of interest, while the lower h curve applies to situations with hori- 
zontal divergence. A later discussion covers procedures used in deter- 
mining which h curves apply for a particular location. The h curve is 
used to determine the mixed-layer thickness in feet from the h scale on 
the left margin of figures 1-8. 

Richard W. James (1966) of the Naval Oceanographic Office has de- 
veloped a new version of thermocline depth prediction based entirely 
on Mazeika's method. This author had completed most of the Pacific 
verification data before James' version was published so that only lim- 
ited testing of it was introduced in this report. 

James believes Mazeika's model to be the most realistic approach 
to wind mixing to date, but points out that no account is taken of the 
demonstrated time lag between onset of wind and the mixing subsequently 
generated, and also that it is overly sensitive to low wave conditions. 
James has corrected these deficiencies in new figures which are described 
below, and has redefined the stability index to what he believes is a 
more realistic parameter. 

James retained the sea state parameter,T)\, as defined by Mazeika, 
but he replaced its determination method by developing a new graph 
(figure 9) which is based on wave forecasting procedures of Pierson, 
Neumann, and James (1955). This can be used in all cases of fully-or 


non-fully-developed sea. 
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To account for the time lag between occurrence of wind and subse- 


quent mixing James developed Table II. This table gives approximate 
time periods for a givenT\ value to stabilize, since a certain wind re- 
quires a given time to develop a maximum7T\ value. If this time require- 
ment is not met then a lower t\ value must be considered. 

James also believed that Mazeika's stability index, AT, is a fair- 
ly conservative factor and that it may not be representative of the re- 
sistance to mixing in many cases. Thus, he redefined the stability 
index (AT, ) as the temperature gradient,in °F/100 feet, at the most 
shallow thermocline. The thermal gradient does not have to extend 


through 100 feet since only the rate of change is desired. 


Salinity corrections for this method are obtained from figure 10. 
The correction is given in °F/100 feet and is added to the previously 
determined stability index (AT,). The amount of correction is propor- 
tional to the salinity gradient, water temperature, and the mixing fac- 
tors. 

James also developed new mixed-layer depth forecasting curves from 


the formula: 


Ki = 
LD = AT, ()— ey Re ATN 


where K and K, are constants. These new curves are depicted in figure 
11. Convergence and divergence was accounted for in this figure by 
using the next AT, curve above the normal AT, curve for convergence, 
and the AT, curve below the proper value for divergence. 

The effects of swell on the mixing of water was also considered 
by James, who states that a much reduced value of T\ should be used for 


swell. The reduction of the mixing force with swell is due to the 
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narrowness of the period band, lack of breaking crests, and reduced 
angular spreading. He also states that swell caused from decaying wind 
waves in a fetch still possesses some of the characteristics of sea and 
should represent about 50 percent of Y\ while swell arriving from a 
distant storm would represent approximately 25 percent of T\. This re- 


duction is strictly a subjective estimate. 
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4.  Mazeika's Prediction Procedure. 

To use the method properly, the stability index must be determined 
first and corrected for salinity gradient if possible. This determina- 
tion can be accomplished in various ways based on: 

(1) averaging temperature differences between the surface and 
the 400-foot level for four to six bathythermograph (BT) observations 
taken within the time interval of 24 hours, preferably; the time inter- 
val could be longer, however, up to a maximum of ten days. 

(2) climatology data records (monthly or seasonal) which give 
sea surface temperature and the temperature at 400 feet. (The accuracy 
of these data will be discussed later.) 

(3) synoptic sea surface chart and climatological data for 
400 feet. 

After AT has been determined, the salinity gradient correction 
can be applied if such data are available. The salinity correction is 
found by entering Table I with the mean temperature in the mixed layer 
and the salinity difference (° Joo) in the mixed layer and coming out 
with a value AT'. This value is then used to determine the correction 
C = - AT'-1° for negative salinity gradients or C =AT'-1° for positive 
salinity gradients, since k(!\) curves actually correspond to 1°F high- 
er stability indexes. This C correction is then applied to the origi- 
nal AT value to determine the proper graph for the mixed layer predic- 
tion. 

The next step in the prediction procedure is to determine the sea 
state parameter T\. This can be done by the use of Table III which lists 


values of wave parameters and minimum fetch and duration for various 
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windspeeds for a fully-developed sea. For application in a non-fully- 


developed sea, H and = are taken from graphs in H.O. 603 and is 


1/3 
then computed with the equation T) = H 


1/3 'max’ 

The forecast windspeed can be a critical factor in determining an 
accurate %\ . The error in mixed-layer depth prediction caused by erro- 
neous windspeed is greater with a small stability index than with a 
larger one. To illustrate, an error of 1 knot in windspeed could pro- 
duce a difference of 16 feet in the predicted mixed-layer depth using 
a stability index of 2° while this same error would produce a differ- 
ence of 6 feet using a stability index of 13°. The accurate forecast 
of windspeed is imperative to an accurate mixed-layer prediction. 

Finally, a decision must be made as to whether an area is under- 
going horizontal divergence or convergence which may influence the ef- 
fectiveness of the mixing processes. This is a difficult problem which 
is approached subjectively through a study of the forecast area. Some 
areas may be under the influence of permanent convergent or divergent 
systems which may preclude conventional mixing processes. To aid in 
this determination Mazeika has illustrated in figures 12 and 13 what 
he believes to be areas of strong convergence and divergence. Figures 
12 and 13 were used as guide lines by this author in the verification 
procedure and they seemed to work quite well. Laevastu and Hubert (1965) 
note that convergence and divergence usually depend on a subjective 
estimation and they recommend the alternative of calculation of diver- 
gence and convergence from surface current fields. 

After the required parameters have been determined and the appro- 
priate graph selected the thermocline depth can be determined. The T\ 


value is found on the far left margin of the graph (figures 1-8). 
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From the horizontal intersection of this value with the proper k(1\) 
solid-line curve (normal, convergent or divergent), proceed vertically 
to intersect the corresponding dashed curve, h(k). Then proceed hori- 
zontally to the h scale on the left margin to determine the thermocline 


depth. 
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a Verification Procedure. 

The requirements for the use of Mazeika's method are explicit. 

To meet these requirements a search for data was initiated and led im- 
mediately to the elimination of all data which included only one BT per 
day. These sources were mainly commercial and military ship reports. 
Since data from Ocean Weather Stations PAPA (50N, 145W) and NOVEMBER 
(30N, 140W) were readily available at the Naval Postgraduate School, a 
decision was made to limit the verification to the Northeast Pacific 
Ocean. Also, the difference in geographical location of these two 
stations in the Northeast Pacific Ocean area would serve to test the 
procedure in perhaps different ocean regions within this broad area. 

Data for Station PAPA were taken from booklets prepared by the 
Canadian Oceanographic Data Center. This source is considered excel- 
lent as all parameters were included in it except the synoptic weather 
picture. Ocean Station NOVEMBER BT data were obtained from the Navy 
Fleet Numerical Weather Facility (FNWF) at Monterey, California. These 
data were in the form of daily reports which normally included from two 
to ten BT observations. 

A point between these two ocean weather stations, which will be 
referred to as MIDPOINT (40N, 140W), was also included. Information 
for MIDPOINT, furnished by FNWF, was available for only one month. 

The synoptic weather charts which were needed to estimate conver- 
gent, divergent or normal synoptic situations were obtained from records 
kept in the meteorology storeroom at the Naval Postgraduate School and 
also from the Synoptic Weather Maps (Daily Series) available through 


the U. S. Government Printing Office. Most of the charts used were 
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facsimile charts sent from the National Meteorological Center at 
Washington, D. C. The particular synoptic chart used included surface 
pressure and 1000-500 MB thickness analysis for the Northern Hemisphere. 
All data collected for the verification were selected from months 
during which the seasonal thermocline is a prominent feature. May 
through August was considered the best period for this to occur at all 
three locations. Station PAPA data included continuous observations 
from May 1964 through October 1964. September and October were included to 
determine the:method's effectiveness after the seasonal thermocline 
period has ended. A limited amount of observations for June 1966 were 
included also. The data periods for Station NOVEMBER included the fol- 
lowing: July 1957, July through September 1958, May and June 1964, June 
1965 and April through June 1966. Data for MIDPOINT were available for 
May 1965. : 
As Mazeika points out in his report, a large amount of continuous 
data (which must include all necessary parameters) from the same observa- 
tion point is not easy to obtain. The data gathered for this verifica- 
tion were found in many different places and brought together so that 
all requirements were met. Numerous other continuous data sources had 
to be discarded as one or two parameters could not be located for the 
time period involved. 
After all available data had been collected, a verification form 
was developed to include all the required information necessary to make 
proper hindcasts. This form is presented in figure 14. Most of the 
column headings are self-explanatory, but explanations for a few are 
required. 
When available, windspeed was obtained from weather observations 
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taken concurrently with the BT observations; otherwise, it was esti- 
mated from the daily synoptic weather charts and entered on the form. 
Much effort was expended in making correct windspeed estimates as large 
errors could develop. The magnitude of this type of error was brought 
out earlier in this report. After windspeed was determined, Table III 
was entered to obtain the value 7%). 

Next, observed mixed-layer depth was entered as the average obtain- 
ed from all BT observations for that day. Since Mazeika's method pre- 
dicts the mean mixed-layer depth it must be compared with an average 
value (since individual BT observations may differ greatly). Tempera- 
tures at the surface and 400 feet were averaged and then compared to 
obtain an accurate stability index. 

After all parameters were obtained, the proper AT graph was used 
in selecting the mixed-layer thickness (thermocline depth). This value 
was then compared with the observed mixed-layer depth. If the observed 
value were larger than the predicted one, the difference between the 
two was entered as a negative number in the Difference column. If the 
predicted value were the larger of the two values, a positive number 
was entered in the Difference column. This was done to determine if 
Mazeika's method was over- or under-predicting the mixed-layer thick- 
ness in the Pacific Ocean. 

When the proper criteria have been met, a check mark is placed in 
the column labelled Verification. In this study a tolerance of +20 
feet from the observed value is considered to have verified. The 
tolerance limits were chosen after considering the factors described in 
the following paragraph. 

First, Mazeika has shown in figure 15 that the average amplitude 
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of internal waves is 17.85 feet. This is interpreted to indicate that 
internal waves alone could provide a variation of thermocline depth of 
this size about its mean value. Second, Tabata and Giovando (1962) 
have shown this same magnitude in figure 16. They also concluded that 
this variation was due to internal waves. Finally, any method which 
can correctly predict within 20 feet of the established thermocline 


depth is certainly a valuable tool for the operating ASW forces. 
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6. Verification Results. 

To insure that Mazeika's method was being applied correctly so 
that the results obtained would be useful, a trial run on Atlantic OWS 
CHARLIE was initiated. A period of one month was tested with the results 
presented in Table IV. The method predicted within 20 feet of the ob-.» 
served mixed-layer depth 63.6% of the time. This was considered accept- 
able, since the hindcasts were made with no previous experience on new 
data, although in the ocean region for which the method was developed. 
After this trial run, it was believed the method was being applied cor- 
rectly. 

Station PAPA had a total of 106 hindcasts performed during the heat- 
ing a May through August. Of these, 82 were verified by the ob- 
served layer depth to show 77.4% of the hindcasts verified within 20 
feet of the observed average layer depth of that day (Table V). In- 
cluding those hindcasts which verified within 30 feet of the observed 
value increased the verification to 85%. Since the mean observed thermo- 
cline depth was approximately 100 feet (Table VI), the above percentages 
can be directly related to verification data presented by Mazeika (1960). 
These percentages are slightly higher than those which Mazeika observed 
on data from Stations BRAVO and DELTA in the Atlantic. However, all 
results obtained from.Station PAPA compared favorably with what Mazeika 
had presented. 

A monthly analysis of Station PAPA verification results presented 
in Tables VII, VIII, and IX shows the high success of the method during 
the heating season. Also shown is a rapid decline in the method's ef- 
fectiveness as the cooling season begins. Table VIII indicates 46% suc- 
cess of prediction was experienced during the month of September and 
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only 21% during October. The period of seasonal or "transient" thermo- 
cline has ended and the mixed layer begins to extend to its permanent 
thermocline depth; this prediction method, which depends on wind and 
mechanical mixing, ceases to be effective since it does not include a 
way of accounting for convective instability in the ocean mixed layer. 

Station PAPA is located in a position similar in some respects to 
the position of Station CHARLIE, at which Mazeika originally developed 
his method. The weather is quite similar in that many pressure systems, 
which normally consist of high winds and associated moderate to severe 
storms, travel eastward through this area. The mixing of the upper lay- 
ers of water due to these disturbances is probably very similar at both 
locations. Also, the latitude of Stations CHARLIE and PAPA is very 
nearly the same. This would indicate that the temperature and radia- 
tion conditions are similar. Acknowledging the fact that water masses 
do vary greatly, the above similarities would lead one to a general ac- 
ceptance that Mazeika's method should be expected to work quite well at 
Station PAPA. 

Returning to Table V, it can be seen that Station NOVEMBER dis- 
plays a completely different set of results than that of Station PAPA. 
NOVEMBER verified only 52 out of 191 attempts for a verification score 
of 27.2%. Table VI displays a mean negative value of error much higher 
than the mean positive value. This indicates that the method was pre- 
dicting a mixed-layer depth much too shallow for the NOVEMBER area. 
Station NOVEMBER is in a semi-convergent mixing zone, so that the area 
is normally, but not continually, subjected to convergent mixing pro- 
cesses. The NOVEMBER area is usually under the influence of high pres- 
sure systems which move eastward. The isobar spacing in these pressure 
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systems is usually quite large resulting in low windspeeds. The low 
windspeed produces a reduced value of sea state parameter which ulti- 
mately leads to a shallower forecast of mixed-layer thickness, using 
Mazeika's method. 

Thus the NOVEMBER area is quite different from the development 
area of Mazeika's method in many respects. The temperature and radia- 
tion conditions differ as well as the windspeed regimes and the overall 
weather patterns. These differences are believed to be the reasons for 
the failure of Mazeika's method in this area. 

With the failure of Mazeika's method clearly evident, it was thought 
a test using James' version would be in order to see if the new stability 
index and forecasting curves might improve the results at OWS NOVEMBER. 
As seen in Table IV the results were also very poor. Only one month's 
data,previously tested using Mazeika's method, were used for this trial 
run. 

Even with the poor results displayed, two interesting side effects 
were noted. First of all, obtaining the temperature gradient of the 
shallowest thermocline did not prove to be an easy task. From the data 
given in various BT reports, the gradient is determined as only a sub- 
jective estimate at best. There were some times, however, when fairly 
accurate gradients could be determined. These occurred when sufficient 
information was reported in the thermocline gradient zone. Secondly, 
the forecast curves of James are a much improved version compared to 
those of Mazeika. They are very easy to interpolate, especially at 
low sea state parameter values. 


Since both Mazeika'’s method and James’ version did not prove to 
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be successful at Station NOVEMBER, one might attempt to improve the 
forecasting there through a new set of curves empirically developed as 
was done at Station CHARLIE. This could also lead to a set of curves 
which would be useful at lower latitudes in general where similar condi- 
tions prevail. 

In the hope of locating a transition area anare the method works 
somewhat better than at NOVEMBER, a point approximately half-way be- 
tween NOVEMBER and PAPA was tested. MIDPOINT had 27 hindcasts perform- 
ed with 8 verifications for a 29.6% predictability figure (Table V). 
This figure indicated that one must be somewhere above 40N latitude in 
order for the method to work more effectively. MIDPOINT also had a 
large negative error value indicating that the method was still not 
forecasting deep enough thermocline depths. 

Tully (1964) describes the behavior of the Northeast Pacific Ocean 
and suggests through his description where the transition between the 
NOVEMBER regime and PAPA regime may occur. In his paper, Tully des- 
cribed various domains or regions which require various models for 
determination of temperature structure and behavior. These domains 
are shown in figure 17. As can be seen MIDPOINT lies in the lower 
part of the Transitional Domain and Station PAPA lies in the center of 
the Central Subarctic Domain. Each domain is considered to have in- 
dividually appropriate meteorological and oceanographic mechanisms 
working within it. If this is true and the study by Tully indicates 
this to be so, then it can be assumed that Mazeika's method should work 
effectively within the Central Subarctic Domain. This is not meant to 
limit the use of Mazeika's method but only to point out an area where 
high predictability reasonably may be expected. 
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In summary, Mazeika's method does produce effective results in 
the Station PAPA area and perhaps these results may be extrapolated to 
include the whole Central Subarctic Domain. At about 45N the method's 
forecast ability may begin to lessen. At MIDPOINT and Station NOVEMBER 
the method is less than 30% reliable, indicating the need for another 
consideration of the forecasting problem in this area; James' version 
of the forecast method did not improve results at NOVEMBER. Mazeika's 
method was easy to apply in most respects except for low values of T\ 
where curves are difficult to interpolate accurately. Some errors in 
forecast mixed-layer thickness may have occurred due to this problem. 

James' version uses a stability index which may be difficult to 
determine in some cases. The thermocline gradient is not always clear- 
ly reported and thus must be estimated. The overall procedure is easy 
to use, especially with regard to James' forecasting curves which offer 


many improvements over the earlier Mazeika version. 
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i= Verification of Climatology Results. 

For the determination of stability index Mazeika suggested as a 
source climatology data (monthly or seasonal) which gives the tempera- 
ture at the surface and 400 feet. Since actual BIT data were used for 
the determination of stability index in this report, a comparison of 
the two was considered for a reliability check of the climatology data. 

Only certain areas have been studied thoroughly enough to pro- 
duce such data. One of these areas is the Gulf of Alaska in the North- 
east Pacific Ocean. Margaret K. Robinson (1957) of Scripps Institution 
of Oceanography conducted a study of sea temperature for this area. 
Robinson analyzed 16,103 individual BT observations from the period 
1941 to 1952. The BT data were subjected to extrapolation, interpola- 
tion and visual smoothing to obtain the various monthly averages. 

Figures 18 and 19 show one version of the temperature averaging 
results. These charts can be used to obtain a monthly average stabil- 
ity index. However, a different version may be more accurate for deter- 
mining temperature from a day-to-day or week-to-week basis. This 
version is displayed in figure 20. This graphical summary was used to 
obtain the temperature data for the comparison. 

Temperatures for the surface and 400 feet were taken from figure 
20 for the months of May through October. These were compared against 
BT observation averages which resulted from this verification report. 
The results have been tabulated and are presented in Table X. 

It must be pointed out that figure 20 is a summary of data for 49N, 
148W which is slightly southwest of Station PAPA. All temperatures 


taken from figure 20 were slightly higher than those which might be 


expected at PAPA, but a comparison should still prove valuable. 
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Since data from Station PAPA were available for the last half of 
May only, Table X includes the averages for this period only. It can 
be seen that the averages compare quite well and in most cases would 
result in using the same forecasting curve for Mazeika's method. These 
averages ought to be even closer if both were done for the same observa- 
tion point. 

From this comparison it is evident that climatology data can be 
used effectively in determining accurate stability indexes in Mazeika's 


thermocline depth forecasts. 
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8. Conclusions and Recommendations. 

Mazeika's method of thermocline depth prediction does work ef- 
fectively at Station PAPA. The method predicted within 20 feet of the 
observed mixed-layer depth 77.4% of the time. This is considered to 
be a good prediction percentage. 

From the study conducted by Tully, it can be assumed the same 
success experienced at Station PAPA may be realized for the entire 
Central Subarctic Domain. This is due to the consistent meteorologi- 
cal and oceanographic properties of this area. 

Station NOVEMBER offered little success with either the Mazeika 
or James version. Both predicted within 20 feet of the thermocline 
depth less than 30% of the time. 

The verification attempt at MIDPOINT resulted in less than 30% 
prediction success also. The success of Mazeika's method at PAPA and 
the failure of it at NOVEMBER and MIDPOINT is primarily due to weather 
and environmental conditions. The temperature, radiation conditions, 
wind regime and weather similarities with regard to Station PAPA and 
the development area of Mazeika's method leads one to believe the method 
should work under these similar conditions. At Station NOVEMBER and 
MIDPOINT the conditions are not the same thus causing the method to 
fail. It is believed that forecasting curves could be developed for 
this area which would make the method useful here. 

During the verification process it was found that Mazeika's fore- 
casting curves were difficult to interpolate at low sea state para- 
meter values. This difficulty may lead to errors in mixed-layer depth 
prediction. All other parameters and procedures were considered very 
easy to use. This method should be adaptable to fleet-wide distribu- 
tion for use in anti-submarine operations. 
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James' version of Mazeika's method was also easy to use, with the 
exception of the stability index. James’ thermocline gradient is more 
difficult to determine than the temperature difference which Mazeika 
uses. The forecasting curves developed by James are much easier to 
use and provide for a more accurate mixed-layer forecast. 

Climatology data which may be used to determine the stability 
index for Mazeika's procedure were also tested. The results indicated 
that climatology data supply accurate information for this determina- 


tion. 
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Figure 20 Graphical summary of BT temperature data at weather station, 49° N, 148° Ww: 
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